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I. HISTORY AND INTRODUCTION 

The U. S. Navy isa hardware oriented organization 
reguiring a wide variety of spare and replacement farts. 
The Naval Supply Systems Command (NAVSUP) is the organiza- 
tion within the Navy responsible for maintaining a sufrfi- 
cient stock of parts and replenishing this stock when 
quantities run low. NAVSUP is keenly interested in the 
quality of the models it uses to provide support for the 
Navy's weapons systems. Faculty of the Naval Postgraduate 
School were asked to examine the existing wholesale and 
retail provisioning and replenishment models and to provide 
comments on their correctness and Suggestions for possible 
improvements. In response, areport by Richards and 
McMasters [Ref. 1] addresses the problems and various solu- 
tion techniyues. This thesis examines a specific solution 
technigue that may be applicable to one or the problems. 
The following is a brief synopsis of the existing models and 
the environment in which they operate. 

Within NAVSUP there are two centers that use provi- 
sioning models, 1) Ships Parts Control Center (SPCC) and 2) 
AVLation /Supply  Cretice a foe Presently the Inventory 
Control Points (ICPs) within the centers have two retail 
provisioning models which generate, 1) COSALS (Coordinateji 
Shipboard Allowance Lists) and 2) AVCALS (Aviation 
Coordinated Allowance Lists). These models supply suffi- 
cient parts support for a short term ship deployment and try 
to provide protection from replacement part depletion until 
the wholesale replerishment models can buy or repair the 
needed spare parts. 

Weapon System procurement requires that provisioning 


support be available by a set preliminary operational 


capability (POC) date. The initial support reguirements are 
determined by a wholesale proviSioning nodel. The support 
must be comprehensive enough to supply a new weapon system 
until the replenishment buy can be made and delivered. si a 
should be noted that the provisioning lead time is usually 
two to two and one half years. Vablaniiley sh product ion 
and installation causes the quality of the forecasted demand 
to be poor. Hence, the proviSioning buy may result in large 
guantities of excess parts in the supply system, or there 
may be many stockouts before the replenishment buy becomes 
available. The Assistant Secretary of Defense 
(Installations and Logistics) has provided guidance [Ref. 2] 
to help prevent large excesses. PROM tn IGm Iinstructzon a 
third model (COSDIF) was developed to conservatively esti- 
mate the quantities of new parts required for wholesale 
provisioning. For existing equipment no additional stock is 
procured for the new weapons systen. If the increase in 
demand from the new system is large enough, the inventory 
manager can anticipate the increase and start a replenish- 
ment buy in sufficient time to meet the need. Otherwise, 
inventory optimization models may be used to control the 
stock. Wholesale replenishment for a new system is initi- 
ated at different times by the two centers within NAVSUP. 
For either center the reorder point for replacement parts is 
subject to change as a function of the actual demand the 
weapoh system places on the supply systen. The reorder 
point for an existing system occurs at part depletion levels 
set by NAVSUP. 

Budgeting constraints for support of existing systems 


are partially driven by the peculiarities of the Department 


of Defense Planning, Programming and Budgeting System 
(PPB). Those purchases for new weapon systems are funded 
as part of the procurement process. For existing systems, 


each Service requests from Congress sufficient funding to 


meet specific levels of readiness and availability set by 
the respective Service Head. NAVSUP receives this direction 
from the Chief of Naval Operations (CNO) andamhen piEvoril 
tizes the parts reguired with a procedure called "Variable 
Threshold". This technique considers several inputs 
peculiar to the weapon system to jevelop the priority list 
for items to be purchased. The quantity (depth) to be 
purchased is determimed from historical data or from engi- 
neering estimates of the number of anticipated failures for 
a period of time egual to one procurement lead time plus cne 
quarter. Proceeding from the top of the priority list the 
parts are purchased to their ‘depth' until the budget 1s 
expended. Those parts unfunded for purchase are placed ina 
file for Unfunded Requirements. 

Okviously, the optimal quantities to purchase, asa 
function of the inputs, and the maximum utilization cf the 
budget are of high interest. Howard in [Ref. 3] addresses 
several different approaches that provide solutions for this 
type problen. The thrust of this thesis is a variation of 
one of those solution technigues for a similarly defined 
provisioning problem. The objective is to minimize the Mean 
Supply Response Time (MSRT) subject to a budget constraint. 

The overall MSRT is a nonlinear composite function of 
the purchase guantity for all items. The individual item's 
MSET functions are also nonlinear with performance as a 
function of four variables peculiar to the item type. ‘The 
budget constraint is a linear function with constant costs 
for each item and a set budget limit. There are no cost 
advantayes for large magnitude buys or item combination 
buys. The individual costs and return functions of the MSRT 
lend themselves to formulation as a dynamic programming 
(DP) problem. However, the problem is large and complex for 
a standard recursive dynamic approach due to the reguired 


huaber of stages, the magnitude of the budget and the range 
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GE costs. consequently, GevarlatiONew OL sw tne Normal D>? 
approach was developed that Slgnificantly reduces’ the 
reguired computations. This variation shall be referred to 
as 'the funnel’ because of its funneling characteristics for 
computation of the resource variable bounds. Further expla- 
NMation of its construction and procedure may be found in 
Clragoe 265 . Modifications of an existing DP computer prograr 
described in [Ref. 3] to operate as the DP variation makes 
the 'funnel't method a usable approach considering speed and 
ease of manipulation. Additionally, the DP variation has 
the advantage of providing an optimal integer solution. In 
this report 'the funnel' will be compared to a full ODP 
procedure and a computerized marginal analysis approach 
which operates very efficiently but which may not generate 
the optimal solution. 
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II. SOLUTION TECHNIQUES 

This chapter will discuss three methods for solving the 
provisioning problem described in Chapter 1, 1) marginal 
analysis, 2) pure dynamic programming, and 3) a modified 
dynamic programming procedure. The first technique provides 
a guick solution that is a good estimate ~but tayeenot be 
optinal. The second technigue provides a global optimal 
solution but reguires a large amount orf computer time for 
any moderate sized problen. The third technigue provides a 
local, and possibly a global, optimal solution and requires 
Significantly less time than option two. 

The following definitions wiil be used in the formula- 
tion and discussions in this chapter. Capital letters 
represent the variable while lower case letters indicate a 
subscript. In the expression Xi-1 the quantity i-1 is the 
subscript although it appears on the same line as the symbol 
Ke 


N = the total number of items considered for 
provisioning 

5 = pre-assigned maximum budget amount 

Cx = the unit cost of item i 

Bae = the essentiality code for item type i 

Ki = the fixed time required to satisfy a 


requisition for item i if stock is 
availakle 
eat = the demand rate for item i(lambda) 
Sa = the decision variable for the number 
of items of type i 
Ld = the procurement lead time for item i 
P'i(Si)= the prokability of Si demands for item i 


during the provisioning intervai 
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Pi (Si) = the probability of Si or fewer demands 
For item i during the provisioning 
interval 

Zi (Si) = the performance measure for item i when 


Si units are stocked. 


The objective of this problem is to minimize the Mean 
Supply Response Time (MSRT) as a function of Si. MSRT is a 
widely used measure cf supply performance because or its 
role in determining availability and because it is an indi- 
cator of the success of a supply system in meeting response 


time goais. It considers two major factors in computaticn, 


1) the likelihocd of satisfying demands from stock on 
hand and 
2) the length of the delay in satisfying demands 


when the system runs out of stock. 


It may be represented by the following expression, 


MSRT=min PEi*Li*Ti*Zi (Si) /PEi*Li*Ti (leat) 
a= Lyre e® @ ai; 


where the likelihood adjustments for conditions one and two 
are factors in the 2Zi(Si) Pune t On. This function repre- 
sents the performance measure of item i when Si items are 
considered and is defined by the following equation for the 
case where demands during the provisioning leadtine are 


assumed to be Poisson distributed: 


Zi (Si) =Kit+((1~Pi (Si)) * (Li*Ti-2Sit Si *(Si+1) 222} 
/LA*¥TLi) + Pt'i (Si) (Li*¥Ti-Si))/2Li). 


The constraint for this problem is the budget or the 
resource in dollars. It is linear and may be represented by 


Eager 2Oltewind tnequality, 
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= Cisa<—e i =teee (228) 


The following subsections present solution techniques which 
are viable procedures for solving the above formulated 


provisioning problem. 


A. MARGINAL ANALYSIS 


The marginal analysis method is a guick, simple and 
effective approach for generating a solution tO the previ- 
ously defined provisicning problem. For the problem consid- 
ered here, the marginal analysis approach generates 
solutions which may not be optimal. However, the solutions 
generated by the procedure do possess some very desirable 
properties which are stated below. 

Let Zi(Si) be the essentiality weighted MSRT for item i 
when Si units are stocked. Let AZi(Si) =Zaatsa- 1) — 2a Si) eee 
the reduction in essentiality weighted MSRT achieved by 
increasing the stockage level for item i from Si-1 units to 
Sie ant Ses Let S'=(S1,S2,...,5N) be the vectoreor i licea— 
tiors for the N items: i.e. Si is the number of units allo- 
cated to item i. The marginal approach begins With a budget 
of B dollars and an allocation vector  S‘=(07 07... 0) Jije 
computes AZi(1)/C1i for i=1,2,2..,N and selects that lteneten 
which this marginal benefit to cost ratio is greatest. 
Without loss of generality , assume that the item selected 
1s iten 1. This reguires an expenditure of C1 dollars and 
results in the allocation vector S'=(1,0,0,...,0). The 
procedure next computes AZ1 (2) /Cl1, the marginal benefit to 
cost ratio achieved by increasing the stockage for item 1 by 
one unit from a level of one to two units. The maximum 
value is then selected from among the ratios in 2.4. 

An additional eunitsot the emcen, say item j, for which the 


maximum occurs is selected and the remaining budget is 
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PVT Omen Bel). se 4 AZNC1) 


Cl ao C3 CN ee 


decremented by Cj dcilars. The process continues in this 
manner allocating a single unit iteratively until the budget 


1s expended or no additional units can be purchased. At 


each step, the ratios compared are 


Le) AZ Se) eee AZN CSN+1). ea 
Cl C2 CN : 


Mieee Sthe Current allocation vector is S$*=(S1,S2,...SN). 
Since only one of the ratios changes at each step of the 
process, the computations required are very few. 

It can be shown [Ref. 4], that for the problem consid- 
epeammnere the solution 1s wndominated. That is, if B{S*") is 
the budget reguired to fund the ailocation vector S', and if 
Seis dM, Other dliocatwon vector such that B(S'*)<=B(S*) 
then the overall essentiality weighted MSRT for S* is 
Blom Ler than Ghat LOLr S* 4. If it should happen that B(S')= 
B, the original budget, for some step of the marginal anal- 
ysis procedure, then the solution S' is optimal. 

In actual practice, what usually happens in the marginal 
analysis procedure is that at some step purchase of the iten 
with the largest marginal benefit to cost ratio reguires an 
expenditure of more money than that which remains. The 
allocation vector at that point is not feasible and various 
heuristics are usually applied to continue the process. The 
most popular heuristic is to backtrack to the previous 
feasible solution and select that item which has the maxinun 
benefit to cost ratio among those items with unit costs no 
larger than the budget remaining. On applying such heuris- 
tics, One Can Not Claim that the final solution so obtained 
is undominated. However, one can make use of the undonmi- 
NMavedsPprEOperty tO Obtdin a bound on the optimal value of the 


objective function. 
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Let Z(SS) be the objective Tfune@menevalicea (ieee a] 
sion at which the marginal analysis procedure first gener- 
ates a sclution for which the budget required is equal to or 
greater than B. If eguality exists, SS is optimal. Suppose 
B(SS)=BS>B and let S*(X) be the optimal solution for a 
budget of X dollars. Then, Since SS is undominated, 
Z(SS) <=Z (S* (BS) ) <=Z2(S*(B)), which 18 tae value oz the ob jece- 
tive function at the optimal solution. Thus Z(SS) isa 
lower bound on the optimal value of the objective function. 
One can thus assess the potential benefit of a more exten- 
sive search for the optimal solution by comparing the actual 
performance obtained at any stajye with the lower bound 
discussed above. 

A computer program was written by Richards to implement 
this method. A sample of this Marginal Analysis prograa 
output in table I shows the iteration number, which items 
ace purchased, the numbers of items purchased, the benefit 


to cost ratios and the budget remaining. The table shows 
the sequence in which the items are purchased. Note the 
decreasing benefit to cost fratio. The program selects iten 


number one five times in a row before the first unit of iten 
two provides a Letter marginal return. The last attempt to 
purchase chooses an item that drives the remaining budget 
negative. The progran automatically re-evaluates the 
Marginal returns and then chooses the next largest benefit 
to cost ratio which in this case drives the remaining budget 
to zero. As previously discussed this solution may not be 
OB tan die 

The Marginal Analysis program has several other capabil- 
ities based on the same approach. There is a job selection 
menu that provides a choice of 1) spares determination, 2) 
budget determination or 3) kit evaluation. The first option 
selects the quantities to be purchased within a given 


budget. The second option solves a related preblem by 
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computing the required budget to meet a specified 
Constraint, 1h this Caseuiek The third option evaluates 
the available performance for a given decision set. The 
program will” conduce @tne Optimization problem with a 
LaGrange Multiplier or marginal analysis. In either case 
the answer is the same. The speed and ease of manipulation 
of this procedure are desirable; however, to be guaranteed 
an optimal solution a means of verification must Fe found. 


That is the goal of this thesis. 


Be STANDARD DYNAMIC PROGRAMMING 


The standard dynamic programming approach can be applied 
to this type of provisioning problen. The following sample 
problem illustrates how the concept may be applied. A 
Single stage problem corresponding to a single item may be 


characterized as in fig 2.1. 





=e ee nn nn 


Y=T(X,S) 





ES ES Te EE Meh te ee 


Figure 2.1 Example of a One-Stage Decision Problen. 


IDerL gure. 2.4, 


X is the input state variable which is the resource 


available for use at this stage. 


S 1s the decision variable that represents the decision 
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within the stage. 


R is the stage return function that produces a return 


based on the input X and the decision S. 
Y is the output state which is a function of X and S. 


Tis the stage transformation function that expresses 
the components of the output state Y as a function 


of the input state and decisions, Y=T(X,S). 


An '*' after an Xi or Si variable indicates the optimal 
Valen ror tidtestage, —~Ref. 5 p2ZZ2=23. }. 


prewone—-Stacge initial-state optimization problem is to find 
the best stage return as a function of the input state X. 
The optimai return from the stage will be denoted as F(X) 


and the cptimal decision as S*(X). Thus, 


F(X)=min &(X,S (X)) =R(X,S* (X)), (2.36) 
S (X) 
where the decision S*(X) provides the best return for that 
particular input value X. If X were to change the values of 
s=(X) and the return function should also change. 

Now consider a multistage problem as shown in figure 
Zz. This is the standard recursive DP structure with back- 
ward numbering where stage Nis on the left and the stage 1 
OnPene right. The optimal return for stages i,i-1,...,1 1S 
represented by Fi(Xi). The solution technique solves 
succeeding stages frem fright to left providing the optinal 
return FN(XN) up to the current stage under consideration. 
In the single stage problem Y was the output state variable, 
but now Y serves as the input to the next stage so 
Moe 1=Y¥=71(Xi;,S1) .- Thus, in an N stage problem, the output 
from stage i becomes the input to stage i-1 for i=1,...,N. 
iieeergure 2.2 


Xi is the state variable 


1g 





Figure 2.2 Example of an N Stage Decision Problen. 


Si is the decision variable 

Ri is the stage return function 

Ti is the stage transformation function 

* after an Xi or Si variable indicates the optimal 


value for that stage, 


where 1=1,...,N is -the stage index. For each stage a deci- 
sion S1i* (Xi) Will be made that yields the optimal return 
Fi (Xi) as da function for the present stage return 
R(Xi,S1* (Xi))} and” the optamal Tcetwrn 1 ee) for the 
Stages, 1 throvgn= 1-4 The recursive eguation for a mini- 


Mization procedure wceuld be, 


Fi(xi)=min (Ri (Xi, Si)+Fi-1(Xi-1)) i=1,....N, (Qe 
BE 
where 

Xi- =T1 CXahg ont) =X1-CiSi Nie eee, Ne (2.8) 


The stage transformation is simply the budget minus the 
quantity purchased multiplied by its cost. This recursion 


assumes the overall return function is the weighted sum of 
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the individual return functions for all replacement parts. 


pmeweguation for this relationship is, 


MS RI=),Ri (Xi, Si) j= 4 eee (2.9) 


The problem is to make a series of decisions SN...S1 that 
will optimize the overall return function while remaining 
within the budget constraints. 

For the previously discussed provisioning problem, using 
the different replacement parts as the stages, the state 
variable Xi is the resource, i.e., the dollars remaining at 
each stage. The decision variable Si is the number of farts 
purchased in stage i. The stage return function is the MSRT 
which must be computed from four operating characteristics 
for each part, (essentiality, procurement lead time, demand 
rate, performance as a function of guantity procured). The 
following equation represents the return function at each 


seage, 


Ri (Xi,Si) =Ei€X¢ Lik Ti*Zi (Si) /_Ei* Lit Ti eon 
1=1,. eee ls 


The Standard DP approach computes the function Fi(Xi) 
meee eCduation 247 Eecursively £0© i=1,...,N and for all 
possible values of Xi. At each stage ithis is done by 
examining all feasible decisions for each Xi, computing the 
return for each decision, and then adding the corresponding 
optimal return Fi(xXi-1) from the remaining stages to arrive 
at the optimal return for the present Xi. Then esOnmall Xa 
values in the stage, the best total return 1s chosen and 
assigned co Fi (Xi). At the Nth stage the process is 
complete and FN (XN) represents the optimal return for the 
problen. This process is efficient for smali problems, 


however it becomes computationaliy more difficult as the 
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number of stages andthe number of possible state variatle 
values increase in each stage. 

For an example of the computational burden, consider a 
very simple problem with a budget 3B, a constant increment K 
for the state variable Xi within each stage, anda max of D 
decisions for each state variable value. In this case tne 
number of evaluations for one stage 1S equal to (B/K) *D. 
Assuning the transformation function and the recursive equa- 
tion for the optimal return equation require only one addi- 
tion each, then there are 2¥*(B/K)*D additions per stage. 
Each state variable value will sreguire D-1 comparisons, 
hence, there are (B/kK) * (D-1) comparisons per stage. The 
total number of operations (additions plus comparisons) 


reyuired for all stages may be represented by, 

N ((2*B*D) /K) +(B¥ (D-1) /K), (a 11) 
or Simplified this becomes 

(N*B/K) *(3D-1) ( 23125 


This example is actually a gross over-sinplification cf the 
computational problen. It 1S possible the return function 
could be very complex reguiring many operations for each 
decision at each Xi value through the range of the budget. 
However, for a numerical comparison, consider a 25-stage 
problem with a budget of 74000, a state variable increment 
of 1000 anda maximum of 100 decisions. This leads to 


953150 operations as shown in equation 2.13, 
( (25* 74000) 71000) * ((3* 100) -1) =553,150. (2. Vey 


This is the required number of additions and comparisons for 


a 25 stage problen. For today's high speed computer an 


oe 


additicn oferation at machine language level takes about 
Beenty—-fLour © machine cycles while a comparison takes 
eighteen. Using twenty cycles as an averaye and a machine 
speed of two million cycles per second, the above number of 
comparisons and additions would take approximately nine 
seconds. This time estimate is for an incredibly simple 


problem that does not take into account any of the adminis- 


trative overhead, ObjeetiveroL Constraint) computations oF 
increased calculations reguired by complex functions. BY 
Setting K=10, which is equivalent to reducing the stage 


increment from 1000 to 10, equation 2.14 shows that the 


humber of oferations will increase to 


((25* 74000) /10) * ( (3*100) - 1) =55,315, 000. (2. 14) 


Hence, it is easy to see the limitation of standard dynamic 
programming. While it gives integer solutions, the computa- 


tional burden is considerable even in simple problems. 


C. A DYNAMIC PROGRAMMING VARIATION 


The computational workload of the full dynamic progran- 
Ming frocedure makes it undesirable for large problems. A 
technigue that will Le referred to as 'FUNNELING! because of 
its control of the resource boundaries between the stages, 
was developed and applied to a standard recursive DP? 
computer progran. The result 1s a shortened procedure that 
provides a locally optimal integer solution to problems that 
can be formulated in the DP construct. 

For any optimal solution to an N-stage DP problem there 
Heya Series of Gecisions Si* £or i=1,;...,N. Corresponding 
to this solution there is a sSeguence of Xi* vaiues for 
lips. sis gitle These optimal state variable values Xi*, repre- 


sent the best quantities of resource availakle for 


ae 


consumption at that stage and are not known in advance. If 


it were possible to find an approximate value for Xi* in 


each stage, the optimal solution then could be localized 
with an iterative prccess. That is, rather than compute 
ECS) for all xX values there could be, a. window "or 


interest, centered around the approximate Xi* value in each 
stage. The computations would be performed on the limited 
range of Xi values in that window thereby reducing the 
required workload significantly. This COnRCeD EeoreconERol 
ling computational boundaries around an approximate optimal 
Xi* value is the foundation for the DP variation. 

The problem initially is how to localize the regions of 
the best Xi. This is solved by the first run of the funnel 
program using a large state variable increment. The 
resulting Xi values are then used aS an approximation to the 
initial Xi* values i=1,...,N. A subroutine of the funnel 
program then constructs a window around the area of interest 
for further refinement. The first estimate of the Xi* 
values 1S guite rough and the decisions, especially the Si 
which have costs less than the initializing increment, may 
change significantly as the Xi are examined with a lesser 
state variable increment. 

The second program run is the first refinement of the 
initial solution. Using the large increments the progran 
will have chosen the best return available though it may not 
be close to the optimal value. The increased sensitivity of 
the reduced state variable increment will allow considera- 
tion of more Xi values in the window permitting the opti- 
Mization process to move closer to the best return. It is 
conceivable that any Xi may change as much as the window of 
its stage will allow. This much variation in the decisicns 
indicates that the best solution lies outside the window and 
the boundaries must be changed accordingly. This situation 


May be present in several stages for any program iteration. 
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Due to the characteristics of the DP procedure, a change in 
the state variable consumed at stage N will cause the 
quantity of resource available at stage N-1 and succeeding 
stages to change. The concept is that repetitive runs of 
the funneling program will allow the computational bounds 
for eacn stage to localize and center upon their respective 
<a When each of the stages has chosen its Xi*, the Si* 
for 1=1,...,N will represent the best solution for a given 
objective and budget. This solution will be a local and 
possibly a global optinuno. Further discussion of the 
program specifics may be found in Chap 3. 

Additional computations may be saved in the funneling 
program by ranking the stages, iargest to smallest, by the 
cost of their respective parts. Purchasing the most expen- 
Sive items first will reduce the magnitude of the budget and 
Simplify the computations in the remaining stages. The 
amount of cemputation depends on three things, 1) the width 
of the window for each stage, 2) the number of state vari- 
able increments within the window and 3) the number of deci- 
Sions available for each state variable increment. 

Consider again the computational example of section 3B, 
where the number of operations is estimated for the dynamic 
programaing solution procedure. To use that analysis for 
the funnel procedure, Bin equation 2.12 is replaced by the 
width of the window W. This width will vary between stages 
but this example will assume Wiis constant and equal to 
10000. Setting the remaining factors the same, K=10, S=100, 
and N=25, the number of operations reguired from equation 


Zeac 1S, 
((25*10000) /10) * ( (2% 100) -1) =7,475,000. (2215) 


This is roughly 13.5% of the full DP requiremeat computed in 


eguation 2.11, a marked improvement in efficiency. There 


22 


are other ways to make the algorithm more efficient and save 
computations. Further discussion of tne details of the 
modified DP procedure may be found in chapter 3. 
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III. THE FUNNEL PROGRAM 
This chapter discusses the operation of the standard DP 
program, the modifications reguired to employ the funneling 
concert and some of the implementation problems encountered 
in the process. Specific subroutine Operations are 
explained in the second section. 


A. MODIFICATION FOR THE FUNNEL PROCEDURE 


The DP variation, called the Funnel program, is a modi- 
fied version of a standard DP computer program named DP5. 
The modified program will be referred to as DP6 for the sake 
of clarity in the following discussion. The DP5 program as 
descrired in [{Ref. 3] has four subroutines, STGPET, STORE, 
tear and DLIM IT. The following list describes what func- 


tion each subroutine ferforms in the DP process. 


1) STORE - enters the constants to be used in the other 


Subroutines and the main progran. 


2) DLIMIT - defines the range of the decision values Si 
for any value of Xi in stage i. 

3) SIGRET - defines and computes the return function for 
any values of Xi and Si in stage i, R1(Xi,51). 

4) TRANFM - defines and computes the transformation 


Pune tncnh tTOrnany Values of X1 and Si in stage 
Providing DP5 with the specific input data as required in 
{fRef. 3], the program will compute the optimal return by the 


standard recursive DP procedure. The following figure shows 


the flow diagram for the general solution technique. 
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Figure 3.1 








Compute recursively 
Fi(X)=min(Ri(X,S)+Fi-1) 
Sa 


TRANSFEM 


| 
using DLIMiIT | 
STGREF ! 


Flow Diagram of the DPS Procedure. 
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The input data are the parameters for the particular problen 
being sclved. When DP5 1S run it uses the STORE subroutine 
to place the parameters in memory, then calls the other 
subroutines, DLIMIT, TRANFM, and STGRET to perform the DP 
procedure. 

To use the concept of funneling described in Chap 2, 
additional subroutines must be added to DP5 to control the 
state variable computational boundaries in each of the 
stages. There are two possible starting conditions for the 
funnel procedure, 1) some input information is available for 
the decisions at each stage, or 2) no prior knowledge of the 
decisions is available. Io respond to these conditions two 
subroutines ENTER and ADJUST have been created to correctly 
Manipulate the inputs. The following list describes what 


function each of the subroutines performs in the Funnei 


program. 

1) ENTER - reads a decision set from a file then 
creates the input data file in the correct 
format for the Funnel progran, 

2) ADJUST - Manipulates the computational boundaries of 


an input data file to reflect the changes 


required at each program iteration. 


The general procedure for DP6 is Shown in figure 3.2, 

The input fparameters are read into memory with the STORE 
fanet lon. Then, depending on the starting condition, the 
ENTER function may used to create an input data file for the 
Funnel program from a previously available solution. Vhs 
solution may be obtained fron, 1) some other solution 
method, 2) a rule-of-thumb estimate of the solution, or from 
3) past information on similar problems. in any of these 
cases the eStimated or best guess decision set is trans- 


formed to the correct input format for the Funnel frogran. 


I: 










is an 

fie. au 
estimate 
avai wabten 











ENTER 
TM taco 
Soult 1 om 











Pertorme recursive, 
Fi (X)=min (Ri (X,S) 
Sa +Fi-1) 
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increment 
and 

windows 
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Baguice 372 Flow diagram of the DP6 procedure. 
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When there is no previous solution, the beginning data file 
may be constructed in the correct format as defined by the 
user instructions located in the comment section of DP6, see 
appendix. Dhis, Ii puree, as described, may be used 
directly with the progran. The next step in the figure is 
the first iteration of the Funnel process utilizing the 
Diet, IRANFU, .and SIGRET subroutines. When the calcula- 
tions are complete the results must be examined to see iff 
the process has found the optimal solution, ies, "LS age 
DONE ‘i. The requirements for the "DONE condition are 1) an 
all integer solution, 2) no window boundaries violated and 
3) the kudget remaining must be less than the smallest cost 
for items to be purchased. If all of these conditions are 
not met the ADJUST subroutine 1S activated and the input 
data file is manipulated to correct the problem areas. With 
each program iteration this corrective procedure is repeated 
until all the conditicns are met and an optimal solution is 
produced. An exampie of the Funnel program output is shown 


fom cigure 3.3. 


| THE OBJECTIVE PSNI 


CONSTBAINT ON COST 


{ 
THE DESCRIPTION WHICH APPEARS BELCW FOR STAGE 1 
APPLIES TC STAGE 1 THRU STAGE 3e INCLUSIVE 
THE PROBLEM IS TO MLINIMIZE A 3 STAGE PROCESSXN 
IS TO BE CHOOSEN OPTIMALLY BETWEEH XN=2.05000D+02 
AND XN=2.0500D+02 
OPTIMAL XN=2-05000D+02 OPTIMAL BETORN=5~30246D+00 
3 2.0500D+02 1.00D+01 5.5000D+01 
2 5.5000Dt01 ~00D+#00 3.50C0D+01 
1 32.5000D+01 7.00D+00 0.0 





Figure 3.3 Example of DP6 Output. 


Bil 


The objective and ccnstraint are stated in the first two 
lines. The description that follows gives the number of 
stages, the budget quantity, the optimal return and detailed 
stage results. Specifically delineated for each of the 
stages are the input budget guantities, the decisions and 
the output budget guantities. For this particular examfle 
the kudget remaining at the. end of the problem is Zero. 

Two ancillary functions were added, one to the main 
program and one to a subroutine to provide informaticn to 
the user. The first 1s a timing function that provides the 
program run time. The second 1S a counting function that 
presents to the screen and stores in a file the numbers of 
the stages in which the optimal .return values are too close 
to the computational boundaries. Both functions have 
provided valuable insights that will be discussed in Chapter 
ae 


Be. THE SUBROUTINE ADJUST 


Ine subroutine ADJUST is the heart of the funneling tech- 
nhigue. To gain the computation reduction desired, this 
subroutine manipulates the windows around the ‘'current' 
State variakle values, Xi', in each stage for each progran 
iteration, Initially the upper and lower limits HXi'" and 
Lxit on the window at stage i are defined by, 


HX1" = Xit+K*Ci+l, (Si. 1) 
and 
LX¥i* = Xi*=K*Gae1, (Be) 


where K 1S a constant. These bounds are illustrated in fig 


3.4. Further adjustment of these limits may be necessary to 


SZ 


assure that one of the Xi values considered in the computa- 
coueone Ftkt) 15 Xa. At stage iin the computation of 


F(Xi) the optimization considers first Xi=LXi, then Xi is 


incremented by Ci to give the next value, etc. fhus 1 2X 
is not properly set, the value Xi=Xi't will never be 
considered. 

pa 

HXi~ | 

| -—— HXi' 

| 

! ! 
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Reece een a! 
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Figure 3.4 Example of a Window on a Stage i. 
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This adjustment of LXi' is made by decrementing from the 
temporary stage value, Xi‘t, by the cost Ci for that stage. 
When this decrementing process first yields a value less 
than or equal to LXi' that value is set as the trnew lower 
limit on the window and is called LXi. A similar adjustment 
is made to the upper limit by incrementing from Xi‘ until 
the value HXi' is equalled or exceeded. The new upper iimit 
is called HXi. This process is repeated for all the stages 
until the funnel is forned. In figure 3.5, the funnel may 
be seen for N stages of a proviSioning problen. The budget 
guantity Bis represented by the vertical lines for all 
stages. The upper and lower bounds are represented by HXi, 
LXi for i=1,...,N where iis the stage index. In each stage 
the bounds are at least a distance equal to the product of 
the constant K and Cit1 from the Xit value. Stages where 
the bounds are furtker from the center indicate the addi- 
tional adjustment for compatability between the successive 


Xi‘ values. 
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The reason for multiplying the succeeding Stage sces a 
the constant K is to allow the DP procedure some computa- 
tional freedom in each of the stages while searching for the 
Xi and Si* (Xi) values. During the computational process, 
should a stage choose an Xi't that is within Cit1 of the 
computational bound, the boundary is said to be violated. 
This is Eecause the DP procedure may actually prefer a value 
outside the bound but is unable to reach it. Figures 3.6 
through figure 3.9 show an example of the stage boundary 
Manipulation as the iterative process of the DP6 program 


searches for the optimal solution. 





—200 =HXi ag i-2 
| 170 =HXi-1 
oe | 120 =HXi-2 
| ol aa | 
~~ oa 
t ==. 20 a2. | 
1100 =Lxi f = ! 
: [70 =LXi-1 Be | 
{ 
20 =LXi-2 
| 
a 20D ea 





Figure 3.6 The Initialization Fun. 


The initializing run depends on the starting conditicn for 
the problem. If the ENTER function is used, the state vari- 
able increments will represent the costs of the respective 
stages. Figure 3.6 shows a starting condition of no prior 
knowledge; hence, the state variable increment is a 
constant, in this example 10, for all stages and of equal 
Humber on either side of the current Xi' estimate. 

A seccnd run of the program provides an opportunity for the 


DP process to evaluate the estimated optimal returns witha 
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Pi gucee3au/ The Second Run. 


smaller state variable increment or a revised boundary value 
PiGdmerOuddjUSt "the windowearound a _omnew Xi' accordingly. 
Figure 3.7 shows the funnel after the second run Eut before 
the windows are adjusted. The program has selected Xi'=140, 
Ma='=70, and Xi-2*=50. These are the current estimates of 
the optimal Xi' values at each stage. Notice that stage 1-1 
PaomamuGindany ss that 1S violated, i.e., the Xi-i' for the 
stage iS within Ci=10 of the boundary. In this case the 
optimization function may prefer an Xi-1 outside the window 
but the boundary prevents its selection. Also in the second 
run the state variable increments have been changed to the 
cost of the item associated with the stage. In most stages 
these values are smaller then the initial increment making 
the return functions more sensitive. This increase in 
sensitivity should improve the optimization selection. 

The violation of the lower boundary in stage 1-1 causes 
the iteraticn to fail the conditions of the DONE checkpoint 
Gumewogune 2.2.) consequently, the ADJUST subroutine is acti- 
vated again to manipulate the boundaries producing the 


results shown in figure 3.8. 
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consumption and avoid cascading boundary violations. 
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Figure 3.9 The Third Run. 
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PeurDe sso shows the funnel after the third run. The 
program has selected Xi'=140, Xi-1'=65 and Xi-2'=46 for the 
current estimates of the Xi‘ values at each stage. Notice 
that the Xi-1i' value is just outside the previous LXi-1, 
which was 70. In this example the Funnel program manipu- 
lated the boundaries of the windows to allow the DP process 
to choose different Xi't values for the three stages, i, i-1, 
ae Ze Subsequent runs of the program will determine the xXi* 
Gon Cach stage. The process is 'DONE' when the remaining 
budget <= mnimum Ci, an optimal integer solution is 


provided, and there are no boundaries violated. 


C. IMPLEMENTATION PROBLEMS 


From the characteristics of the funneling process, there 
are prebliems that arise when transitioning between stages. 
When computing the maximum and minimum value of the decision 
to be made for each Xi in a stage i, the upper and lower 
bounds and item cost cf stage 1-1 must be considered. 

In figure 3.10, the left side is the window for stage i, the 
Pitiat siae 1S the window for stage i-1. For the stages i, 
i-1, the highest values of the state variables are repre- 
sented by HXi and HXi-1 while the lowest values are repre- 
sented by LXi, LXi-1. In the program DP6, the two variables 
DLOW and DHIGH represent the mininum and maximum decisions 
availabie for the state variable in each stage. eS 
essential that DHIGH for any Xi value in the window of the 
Htmstage be Sf such a Magnitude that it will not force the 
Xi-1 value relow the lower bound LxXi-1. AGG@lironay.y, DLew 
must be of such magnitude that for any Xi in the window of 
Giemremestade, it will not. Eorce the Xi=?o value to violate 
the upper bound HXi-1. If either of these decision quanti- 
ties is incorrect at stage i the resource quantity will be 


outside the window at stage i-1. The derivation of the 
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Figure 3.10 Boundary Manipulation. 


formulas for these values follows. The transitormatton Une 


tion as represented by equation 2.8 of Chap 2 is, 


Xi-1=Xi-CiSi i=1,...,N. meh 


If this transformation is to take place from a window in 
stage i to awindow in stage i~1 the following conditions 


must be true. 


Di 1 <= - Cre Ga, (3.4) 


where SHIGH is the maximum number of items to be purchased 


and 

Xi-Ci*SLOW <= Hxi=—1, (3.5) 
where SLCW is the minimum number of parts to be purchased. 
By algebraic manipulation a relationship for SLOW and SHIGH 


may be established, 


4Q 


SHIGH <= (Xi-LXi-1)/Ci, (3.6) 


and 


SLOW >= (Xi-HXi-1) /Ci. (Ss) 


The expressions SHIGH AND SLOW are equivalent to the DHIGH, 
DLOW required in DP6 and are coded in the subroutine DLIMIT 
to control the boundary calculations. It is conceivable 
that at the extreme state variable values of Xi the computed 
boundaries could become negative or larger than the assigned 
budget. Thus, the decisons are further limited Ly DHIGH= 
min(SHIGH,200) and DICW= max(0, SLOW), where 200 is an arbi- 
trarily chosen value, 

The initialization for the problem is either from input 
decisicns or else a large state variaple increment is used. 
After the subroutine ADJUST manipulates the windows the 
first time, the state variable increments are decreased or 
the boundary conditions are changed to allow examination of 
the new possible Xi values for the best return. The subrou- 
tine DLIMIT automatically adjusts the DHIGH and DLOW values 
to compensate for the changes in the state variatle incre- 
ment. With this closer scrutiny, the optimal return value 
may change. The ADJUST subroutine moves the boundaries up 
or down accordingly when an optimal value is within one 
ancremeat of the bounds. Recall that if the optimal return 
value is within one increment of window bound, it may be 
that the DP process would prefer a value outside the window. 
Hence, the boundary is recorded as violated, and the window 
is moved for the next iteration. This process is repeated 
until all the optimal return values are centered within the 
windows for all stages. 

TO save computations and avoid computing Xi values that 


are not compatible between stages, the state variable 
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increment for each stage is the cost of the item to be 
purchased in the previous stage. A problem arises in that 
the original DP5 program may frovide solutions tnat are not 
integer. This is fossible because in the stage being 
considered, the best return may reyuire an Xi't value fron 
succeeding stages that has not been computed. DP5 handles 
this by interpolating linearly between the state variable 
values it does have to arrive at an estimate for the one 
required. In DP6 tc avoid thiS approximation and oktain 
integer solutions, the subroutine ADJUST, (as discussed in 
Section A), manipulates the lower and upper window values to 
ensure the previous stage optimal return values are avail- 
able for the Xi being considered. This adjustment process 
is also repetitive and is done in conjunction with the 
boundary control. 

Another problem arises from the effort to save computer 
computation time. If there is a wide range of costs, the 
use of the item cost as the state variabie increment within 
a stage becomes troublesome. Initially the ADJ US? s wero 
tine creates the state variable windows for the individual 
Stages using a constant multiplier that may be set by the 
user. To remedy the problem of the extreme range of the 
-costs of items for the proviSioning problem being consid- 
ered, ($10000.00 to .50), the multiplier was changed to a 
variable that depends directly on the item cost. This was 
done to decrease the program repetitions reguired to reach 
an oOptinal scolution. Under the constant multiplier method, 
a change inthe purchase guantity of a larger cost iten, 
(i1.e.$1500), could reyuire many program repetitions to 
re-allocate the resource available for the lower cost items, 
(1,e.310,5,3.75, etc). The smaller window range, gue stceurn— 
constant multiplier, hampers the re-allocation effort. The 
variable multiplier allows greater compensation for lower 


cost items to facilitate a more rapid change if required. 
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IV. ANALYSIS 

The primary topic of this chapter will be the comparison 
or the Marginal Analysis and Funneling solution technigues 
for the previously stated problen. The discussion will 
cover the results of both programs and the solution differ- 
ence. The standard dynamic program procedure will not be 
addressed because orf the computational advantage of the DP 
variant and the similiarity of the two techniques. 

For Simplicity and ease of comparison a small _ three 
stage problem will be defined and solved using both solution 
procedures. As in the original provisioning problem the 
opjective is to Minimize the MSRT with a constraint on cost. 
Three different budget guantities will be examined to denmon- 
strate the differences in the solution technigques. The 
following table is a list of the input parameters which are 
cost Ci, demand rate Li, essentiality El, and procurement 


lead time Ti. For this particular example Ti is one year. 


TABLE If 
Sample Problem Input Parameters 


i A yD ce ee ee ee ee 


| 

| I DEMAND COST T ESSEN ) 

| 1 5.000 5.00 1.00 3.0 | 
2 1.000 10.00 1.00 oO | 
3 10.000 16.00 1.00 1.0 | 

ee ee ee amie a eed 
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For the first budget quantity of $195.00 the Marggnal 


Analysis program returns the following output. 


ea SS ON ae MS a EE 


TABLE IIifi 
Marginal Analysis Program Return, Budget=$195.00 
EHEKKEEKHEKE RAKE EEA EREA EEE EAEE REEDS SEES E ET 
+ METHOD Ser ear BUDGET = 195.00 * 
EKKEKTKKEKREKSEKKKEKSEKAKEKRKKSE ESSERE AEE RKKKHKEKEFELESES 


PCC CTS ECE SCESE SESES SE SEES ESTES ESSERE SES TSE SESE SE EE 


STEP ILEe ot BATIO BR 
i 1 1 0.480808496 190.00 
2 1 2 0.365659833 185.00 
3 1 3 0.260€18091 180.00 
a 1 4h 206 17 2527 175.00 
5 1 5 0.105280340 170.00 
6 ie 1 0.0735758543 160.00 
q 3 1 0.0600002967 145.00 
8 1 6 O.05979566/5 140.00 
9 a 2. 0.053 35695 93 125.900 
t0 3 3 0.0466887541 110.00 
11 3 G4 0.0400909968 95.00 
i2 3 5 0.20336193480 80.00 
13 i 7 0.030657734) 75.00 
a4 3 6 0.0273999237 60.00 
15 3 7 0.02160086488 45.00 
16 2 2 0.0207276568 35.00 
17 3 BS 020 160023277 20.00 
18 q 8 0.0146531016 15.00 
19 5 9 0.0131195445677 0.0 


RHEE EK ERK RE ERERARAER ERE AEE RETHEEEERES 
MINIMIZE MSET EXCESS: 0.90 


(ines: cee Mn ees ace ci ea me (A mia _— 


I Si MSRT , 

1 3 1,302 ! 

2 2 10.396 ! 

3 9 122196 | 

| 

OVERALL 6.003 | 
REALIZED PERFORMANCE = 6.0033 | 
BOUNDS { 6.0033, 6.0033) | 
| 


SS SS ee cs SS OS 
ee ee ee ee ee 


The decision results are, $1=8, S2=2, S$3=9 and the return is 
Z= 6.0033, with no kudget remaining. Viewing the purchase 
sequence intable III it may be seen that the Marginal 
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Aualysis prcgram has reached a solution uninterrupted, that 
is without implementing the heuristic developed in [Ref. 1]. 
Consequently, for this budget quantity tne solution is 
optimal. Notice ancther indicator of the optimal solution 
in the last line of table III. The bounds described are the 
possible upper and lower values for the return. When these 
Guantities are the same value, the solution is the best 
available. 


: 
| 


| | 
TABLE IV | 
Funnel Program Results, Budget=$195.00 
THE OBJEGZIVE IS SSR | 
CONSTRBAIST ON COST | 
| THE DESCRIPTION WHICH APPEARS BELOW FOR STAGE 1 | 
APPLIES TO STAGE 1 THRU STAGE 3 INCLUSIVE 
THE PROBLEM IS TO MINIMIZE A 3 STAGE PROCESSXN | 
IS TQ BE Been. OPTIMALLY BETWEEN XN=1.9500D+02 
AND XN=1.9500D+0 
OPTIMAL XN=1.95000D#02 OPTIMAL RETURN=6.00332D+00 | 
4 XN S XN—- 1 
3 1.9500D+02 9. 00D+00 6.0000D+01 | 
2 6.0000D401 2. 00D+00 4. O000D+01 | 
4.Q0000D G20 | 
mao 


e 
| 1 +01 8. 00D+00 





The same budget and parameters entered in the Funnel 
program produce the results in table IV. The output is in 
the same format discussed in Chapter 3. The initial budget 
quantity, XN and the optimal return are shown in the line 
just above the individual stage results. The variable Xi is 
the budget guantity at the beginning of the stage, S1 1s the 
decision or quantity purchased, and Xi-1 is the budget at 
the end cf the stage. Recall that item 3 is the most expen- 
Sive and item 1 the least expensive. The optimal solution 
for the $195.00 budget is $1*=8, S2*=2, and S3*=9 with a Z 
value equal to 6.0033 and no budget remaining. This answer 
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matches exactly the Marginal 


Analysis program results 


shown 


in table III, confirming that the Marginal Analysis solution 


IS¥Opti was 


- 


a et a. 


Was 


Marginal Analysis Program Return, 
ee ARH AKAEKEEKEREAERARAAERE KAAS SIGS EKER RED 
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ETHOD: FIXxt 
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ABLE V 


Budget=$200.00 


UDGE 200.0 
EREKHKEEKAKHKREEKHKAAHK RA EKARERERE RAAT ERAKEE THEE DF 


—o 


* 
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pS SPESSEESISOLSESSLSSSESSESSEESE LSE SESE ESSE SEE SE SS SE ® | 


MINIMIZE MSRT 


R 
B 


EA 
OU 


1 
: 
2 
3 


Li 
NDS 


FOr the 


used. 


SI 
9 
2 1 
9 1 

OVERALL 


ZED PERFORM 
( 6.003 
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EXCESS: 
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budget quantity 


eet cen St pment ates Sc ee ena aan anes SO a GANED cnet EnDhenanms ste gps, SSID natn eres si PU cena PU ga fNS acta mut corre AOID gene eth grape ca ge ms SO eden Gee see 


of $200.00 


The Marginai Analysis program returns the output 


in table V. The results are S3=9, S2=2, and $S1=9 with a Z 
Maree Ge 5.5052 and no budget ieft over. Examining the 
detailed output of the Marginal Analysis program, table V, 
it may be seen that the purchase of the twentieth iten 
Grives the budget negative. The program impiements the 
heuristic, re-evaluating the benefit to cost ratios, then 
chooses the humber one item for the next purchase. Notice 
that in the bottom line of tabie V the upper and Jlower 
bounds are not equal. [vcs weomeamgooce indication the solu- 


tion may not be optimal and should be verified. 


ee a Sl a i ge a 
TABLE VI | 

Funnel Program Results, Budget=$200.00 
THE OBJECTIVE IS MSRT | 
CONSTRAINT ON COST . 
THE DESCRIPTION WHICH APPEARS BELCW FOR STAGE 1 
APPLIES TO STAGE 1 THRU STAGE 3 INCLUSIVE ! 
THE PROBLEM IS TO MINIMNIZE A SVSLAGE _PROCESSXN | 
i IS TO BE CHOOSEN OPTIMALLY BETWEEN XN=2.00C0D#t02 j 
AND XN=2.0000D+ 
OPTIMAL XN=2.0000D402 ene OE eo ee } 
N XN i~ : 
: 3 2.0000D402 9.00D+00 6.500CD+01 | 
2 6.5000D+01 2.00D+C0 4.5000D+01 
1 4.5000D+01 9.00D+00 0.0 { 
— eee eee eee es: 


The same budget and parameters entered in the Funnel 
program produce the results yiven in table VI. In this case 
the answer is again identical, $1=9, S2=2, and S3=9 with the 
optimal return Z= 5.56519. Incidentally, the heuristic in 
the Marginal Analysis program makes the correct choice at 
the twentieth step when it selects the item with the second 
highest benefit to ccst ratio. 

PoOmeameiind Case, dunudges: of $205.00 1s entered in Doth 
programs. The resuits for the Marginal Analysis code are 
Shown in table VII. 
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TABLE Wil 


Marginal Analysis Program Results, Budget=$205.00 
KEKEKEAKEA AEE EERE EERE EERE AERA EKA TEREE REEDED 
* MET HCD: FILED EUDGET = 205200 * 
REKKEKEREAEREKEKERER EERE RARE REAR EDER OKT EEE EEE ES 


KKKKKKKKEKHKKKKEKERSEKSEEKSEKEKKAEKESHKKEKEEKRESLERESRS 


© 
ty 
= 
3 
— 
© 


STEP ITEM SI B 
1 1 1 02480808496 200.00 
2 1 2 0.365659833 195.00 
3 1 3 02260618091 190-00 
4 1 G 0. 172421157 185-00 
5 1 5 0.105280340 180.00 
6 2 1 0.0735758543 170.00 
7 3 1 0.0600002967 155.00 
8 1 6 0.0591956675 150.00 
9 3 2 0.0533369593 135.00 
10 3 3 020466837541 120.00 
ile 3 4 0-0460909968 105.00 
12 3 5 0.0336193430 90.00 
13 1 7 0.0306577347 85.00 | 
14 3 6 0.0273999237 70-00 
15 3 7 0.0216008648 55.00 | 
16 2 2 020207276568 45-00 | 
17 3 8 0.0164023377 30-00 
18 1 8 0.0146531016 25.00 
19 3 9 0.0119544677 10-00 
20 3 10 0.0083406679 -5.00 
20 1 9 0.00648 18673 5.00 
21 1 10 020026624922 0-0 


KEKE EERE EEK RASH HE REE REET EEEEAEK EDS 
MINIMIZE MSRT EXCESS: 0.0 


I oe Sma 
1 10 0. 190 
2 ia 10.396 
3 9 12.176 
OVERALL 5.385 

AN 

34 


ood (a) ened 


‘men fa) fag ee ae Pcie) amend Feira a fe a I im deme ie eG enable 
=_—=a emi pci ies Elite 


REALIZED PERFORM a6 a2 


= Dis 
BOUNDS { 6.003 4.3120) 


aha Dever onan 5 Slants eadbad ah 9 Ba a 


The solution is S1=10, S2=2, and S$3=9 with a return value of 
i= Seduce The upper and lower bounds are not equal indi- 
cating the heuristic has been applied and the solution 
Should be verified. Examining the output from the Marginal 
Analysis program, tablesvpi, at May be seen the budget was 
again driven negative at the twentieth item. The heuristic 
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waS activated purchasing two of the number one items to 
drive the budget to zero. 

Entering the same input and budget parameters in the 
Funnel program produces the results shown in table VIII. 


Funnel Program Results, Budget=$205.00 


— -— at ema eg A A A RR RS ae FG GA PSR RR RE AS SE RANA SERRANO at MEIN 


TABLE VIII | 
| 


THE OBJECTAVE IS MSR 
CONSTRAINT ON COST 


= = 
a 
| 
| 








THE DESCRIPTION WHICH APPEARS BELOW FOR STAGE 1 
APPLIES TC STAGE 1 THRU STAGE 5 LNCLUSIVE 
THE PROBLEM IS TO MINIMNIZE A 3 p51 AGE PROCESSZN 
IS TO BE CHOOSEN OPTIMALLY BETREER XN=2.05000D+02 
AND XN=2.0500D+02 | 
es oe Seen Ene uO ot OC 
2.0500D+02 1.00D+01 5.5000D+01 | 
5-5000D+01 - 00D+00 ~20C00D+t01 | 
325000D¢01 - 00D+00 0.0 
The answers are not the same. The Funnel program resuits 
indicate the best part selection is Sl=7, S2=2, and 5S3=10 
for an optimal return of Z= 5.30246. This selectior 


provides an MSRT that is roughly .08 better than the 
marginal analysis method. 

Recall from section A of Chapter 2 that the Marginal 
Analysis program has a function called KIT that will deter- 
Mine the return froma given decision set. Entering the 
optimal solution from the Funnel program into the KIT func- 
tion produces the result in table IX. The resulting objec- 
tive value is identical to the value computed by the DP 
Variant. This shows that the computational procedures are 
the same for both programs. 

Now consider the definition of a much larger problen. 


It is of the same type as the previous examples, an 
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TABLE IX 
Kit Function Results, Budget=$205.00 


EVALUATION OF SPECIFIED ALLCCATION.... 


1 SI MSRT 
1 7 3.085 
2 iz 10.396 
3 10 72610 
OVERALL Se 2 
REKKKAK KEKE EKER HERA EESR ARERR DEE E EEE EERE ESE D 
~ ee 





objective of Minimizing MSRT with a constraimeyon Gost nc 
input parameters cost Ci, demand rate Li, essentiality fi, 
and procurement lead time Ti, are shown in in table X. The 
budget is $74825.00 with twenty five stages and costs 
ranging from 3170000200 “tomo Due to the length of the 
detailed results an abbreviated output of the Marginal 
Analysis solution is shown in table XI. The return is 2. 853 
With this decision set. Note in the last line of the table, 
the split in the bounds indicates the heuristic was imple- 
mented and the solution may not be optimal. Examining the 
abbreviated output of the Marginal Analysis program, MTakle 
3.4 shows that the budget was driven negative ten iterations 
from the end. Entering the same parameters and budget gquan- 
tity for bcth starting conditions of (the Funnel S@peeqran 
produces identical results shown in table (XII. Comparing 
the results in tables XI and XII it may be seen that the 
return and decision set for the problem are the same for 
both programs. Interestingly, the Marginai Analysis program 
was able to correctly select ten items after the heuristic 
implementation to match the optimal solution produced by the 
Funnel fprogran. Tke heuristic employed by Richards and 
NcMasters appears to have been a good choice, although the 
solution must still be verified. 
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TABLE X 
Input Parameters for a Larger Problem 
COST 
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BABLE AL 


Marginal Analysis Program Results, Budget 


$74825.00 
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V. CONCLUSIONS 

This chapter will discuss the advantages and disadvan- 
tages of the two primary solution technigues examined in 
this paper. Various starting modes of the Funnel program 
will be looked at in an effort to develop an efficient 
approach to this type of provisioning problen. 

The Marginal analysis program provides a guick, simple 
and effective approach for jenerating a solution to the 
previously defined provisioning problem. However, the solu- 
tion is not guaranteed to be optimal because the prcgram mav 
implement the Richards heuristic to consume ail avaiiatle 
resources. The Funnel program provides the optimal solution 
to the same problem though it takes more computer time to 
operate. The DP variant has the capability to start the 
problem from two different conditions, 1) no prior knowledge 
and 2) prior knowledge or an estimate of the solution. For 
Case 2, utilizing the ENTER function of the Funnel progran 
allows the rapid construction of the input data file used by 
DP6 to examine the Xi* values and manipulate the window 
boundaries as necessary. If the entered solution set is not 
optimal, the Funnel program will drive the decisions to the 
optimal decision set in subseguent iterations. Using the 
timing function mentioned in Chap. 3 a rough approximation 
of run time is available for each iteration of the progran. 
For Case 1, no prior knowledge, the total estimated time for 
the $74825.00 problem is nearly 3000 seconds. FOr CaseuZ, 
prior knewledge, a near optimal starting solution for the 
same problem is driven to the optimal solution in apfroxi- 
mately 300 seconds. This 90% savings on computation tine 
Makes it clear that an approximation method as an input to 


the Funnel frogram would be a better approach to this type 
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Ole roblen. It should be noted that the run time for the 
case 2 start is related to the accuracy of the entered deci- 
sion set. In this instance the entered decision set was 
Specifically constructed to force the Funnel program to 
adjust the window boundaries in all stages at least once. A 
better solution than the one entered could produce a muck 
shorter run time. 

There are tradeoffs in the Funnel Program for computa- 
tional efficiency. Variations of the multiplier that 
creates the window can drive the individual program run time 
directly. A window that is too larye wastes computation 
time on resource values that are not needed. tf the window 
width is too small, multiple iterations of the program may 
be reguired to compensate for the increased number of 
boundary violations in the effort to reallocate the resource 
consumed at the various stages. HORmEne Verliitcation of a 
deciSsicn set provided by the Marginal Analysis program the 
ioe ner COUuLd be set to 1.0 providing avery narrow 
window of resource to be examined. Due to the narrow width 
of the window, if the Xi* value changes in any stage the 
boundary will be recorded and displayed to the user as 
violated. This notification would indicate that the entered 
solution is not optimal for the problen. 

The computer time expended for the verification alterna- 
tive would te quite low even for large budget problems. In 
the 25-stage problem, for example, approximately one second 
is reguired to verify the correctness of the optimal solu- 
Biome etters Entered initially. if the verification alter- 
native indicates the decision is not optimal, the time 
reguired to drive the entered solution to the optimal solu- 
tion is controlled by the magnitude of the changes in the 
decision set and the number of iterations reguired to 
re-allocate the resources throughout all the stages. The 
user must make a value judgement on the importance of the 


optimal solution versus the additional computer costs. 


Sle, 


In cenclusion, the emphasis of this paper 1S to -@beate 
an algorithm to determine the optimal solution in a quick 
and easy manner. It is evident that a combination of the 
Marginal Analysis and Funnel programs is an efficient 
approach for solving this type of problen. Using the 
Macginai Analysis program to generate a best estimate of the 
solution, the ENTER function with the Funnel program can 
then ke applied to ccnfirm the optimal solution or to create 
an input data file that DP6 may use to find the optimal 
solution. The provisioning problems faced by NAVSUP, as 
described in Chap 1, may be solved with this method if they 
can be formatted in the recursive DP construct. The prograa 
DP6 in particular solves only the minimization of MSRT with 
a constraint on cosi. Variations on this program could be 
devised to solve additional objectives like supply material 
availability (SMA), time weighted units short(TWUS), etc., by 


Changing the computation procedures in the subroutines. The 


constraint however, 1s currently restricted to the linear 
Toren. 
EONIAS1<=Baee— 1,25 20 (5.1) 


where Wi is resources required per unit, Si is the number of 


units of item i, and B is the total resource available. 


a0 


SS SoS eee SS Se 


DP6 CODE 


This appendix contains the computer listing of the main 
program and subroutines for the funneling procedure. Some 
of the Subroutines contain code for objectives and 


constraints other than those discussed in this thesis. 
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